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Abstract. An operator H = Hq + V where Hq = i~^d^ {N is arbitrary) and 
y is a differential operator of order — 1 with coefficients decaying sufficiently 
rapidly at infinity is considered in the space L'^(IR). The goal of the paper is 
to find an expression for the trace of the difference of the resolvents {H — z)~^ 
and {Ho — z)~^ in terms of the Wronskian of appropriate solutions to the 
differential equation Hu = zu. This also leads to a representation for the 
perturbation determinant of the pair Ho,H. 



1. Introduction 

1.1. In the framework of the general operator theory in an abstract Hilbert 
space, the spectral theory of differential operators 

H = i-'^d^ + VNix)d'^-^ + ■■■ + V2ix)d + viix), d^d/dx, (1.1) 

is the same for all values of N. However, from the point of view of differential 
equations the problems are essentially different for = 2 (for A'' = 1 it is trivial) 
and for larger values of N. 

Suppose that the coefficients Vj{x), j = I, . . . ,N, decay sufficiently rapidly as 
|a;| ^ oo, and set Hq = i'^d'^. Let Rq{z) = {Hq - z)-^ and R{z) ^ {H - z)-^ 
be the resolvents of the operators Hq and H acting in the space L^(M). The self- 
adjointness of the operator H is inessential for us, and we do not assume it. 

The main goal of the present paper is to find an expression for the trace 

Tt (Riz) - Roiz)) (1.2) 

in terms of solutions to the differential equation Hu — zu. In the case N = 2 such 
an expression was found by V. S. Buslaev and L. D. Faddeev in paper [S]. They 
considered the problem on the half-line, and the problem on the whole line was 
discussed by L. D. Faddeev in [7]. 
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1.2. Let us introduce the notation 

I Ul ... UN \^ 

{ui,...,ujv} = 



(1.3) 



for the Wronskian matrix of solutions ui = ui(x, z), . . . ,un — un{x,z) of the 
differential equation 

i-^u'^^\x) + vn{x)u'^^-^\x) + ■■■ + V2{x)u\x) + vi{x)u{x) = zu{x). (1.4) 

We always assume that 2: G C \ [0, 00) if N is even and that Imz ^ if TV is odd. 
Let C,j be solutions of the equation — z. We suppose that 

ReCj>0 for j = l,...,n and ReCj<0 for j = n + 1, . . . , iV. (1.5) 

Here n = N/2 if N is even and n = (iV - l)/2 for Im z > and n = (A^ + l)/2 for 
Imz < if TV is odd. 

We first explain our result for the case of functions Vj{x) with compact supports. 
We write cc « if x lies to the left of the supports of all Vj (x) and x :» if x lies 
to the right of this set. Let Wj(x, z) be solutions of equation (|1.4I) such that 

Uj (x, z) — e'^'^ for X <K if j = 1, . . . , n and for x » if = n + 1, . . . , A^. 

(1.6) 

Let 

W(x, z) = det{ui(x, z), . . . , u„{x, z), -u„+i(x, z), . . . , un{x, z)} (1.7) 
be the determinant of matrix (II. 3L and let 



Wo(z) = det{e^i", . . . , e«"", 6^"+^", . . . , e^""} (1.8) 

be the corresponding Wronskian for the "free" case where Vj — for all j = 
1, . . . ,N. Of course, 

W(x2, z) = exp ( - /" ' VN{y)dy)\N{xi,z) (1.9) 

for arbitrary points xi and X2. We emphasize that the Wronskians \N{x,z) and 
Wo(z) depend on the order of numeration of the numbers but the normalized 
Wronskian 

A(x, z) = W(x, z)/Wo(z) (1.10) 

does not depend on it. 

Our main result is that the normalized Wronskian satisfies (for all x and all 
regular points z of the operator H) the equation 

Tr (i?(z) - i?o(z)) = -A(x, z)~^dA(x, z)/dz, (1.11) 

which we call the trace formula in this paper. Thus the trace of the difference of 
the resolvents admits an explicit expression in terms of properly chosen solutions 
of equation (|1.4p . 

Then we extend representation (|l.lip to general short-range coefficients Vj{x) 
satisfying the assumption 

|wj(x)p(l + x2)"dx < 00, a>l/2, j = l,...,N, (1.12) 
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only. In this case the functions Uj{x,z) in definition (|1.7p are the solutions of 
equation (|1.4|) such that 

Uj{x,z) = e^'''{l + o{l)) (1.13) 

as X — oo if J = 1, . . . , n and as x — +cxd if j = n + 1, . . . , A^. Here and below 
all asymptotic relations for solutions of equation (|1.4p are supposed to be N — 1 
times differentiable in x. We emphasize that for N > 2 asymptotics (|1.13|) DO 
NOT determine the solutions of equation (|1.4|) uniquely. However, the Wronskian 
(|1.7p does not depend on specific choice of the solutions satisfying (|1.13p . Thus we 
do not need the construction of the book [2] by R. Beals, P. Deift and C. Tomei 
devoted to the inverse scattering problem. In [2| solutions of equation ()1.4p were 
distinguished uniquely (away from some exceptional set of values of z) by conditions 
at both infinities. Our construction of solutions of equation (|1.4|) with asymptotics 
(jl.l3p relies on integral equations which are Volterra equations for N = 2 but are 
only Fredholm equations in the general case. Nevertheless for the construction of 
solutions with asymptotics (|1.13l) as x — > +cxd (as x —oo) we impose conditions 
on the coefficients Vj{x) also as a; — ;> +oo (as x — — oo) only. 

Suppose that vn — 0. Then W(a;, z) = W(z) and hence A(x, z) — A{z) do 
not depend on x. In this case we identify A(z) with the perturbation determinant 
for the pair of operators Hq, H. We refer to the book [TO] by I. C. Gohberg and 
M. G. Krem for a comprehensive discussion of different properties of perturbation 
determinants. Set 

V = H-Ho^ WAr(x)a^-i + • ■ • + V2{x)d + vi{x). (1.14) 

If Vn — 0, then the operator VRo{z) for Imz ^ belongs to the trace class 6i, 
and hence the perturbation determinant 

D{z) = Bet (I + VRo{z)) (1.15) 

is well defined. Of particular importance is the abstract trace formula 

Tr {R{z) - Ro{z)) = -D{z)-^dD{z)/dz, (1.16) 

which for definition p.lSp is a direct consequence of the formula for the derivative 
of a determinant. Comparing equations (|l.lip and (|1.16p and using that A(z) 1 
as I Imz| oo, we show that 

Det (/ + VRo{z)) = A(z). (1.17) 

Thus the perturbation determinant admits an explicit expression in terms of solu- 
tions of equation (|1.4p . 

li Vn 0, then under assumption (|1.12l) it is still true that (for all regular points 

^) 

Riz) ~ Ro{z) e 6i, (1.18) 
although VRo{z) ^ ©i. Without the condition vn = 0, equation ()1.16p is satisfied 
for so called generalized perturbation determinants D{z) which are defined up to 
constant factors (see subs. 6.2). According to equation (II. lip in the general case for 
every fixed x G R, the function A(x, z) differs from each generalized perturbation 
determinant by a constant (not depending on z) factor. 

1.3. A preliminary step in the proof of the trace formula (jl.lip is to find a con- 
venient representation for the resolvent R{z) of the operator H. This construction 
goes probably back to the beginning of the twentieth century. We refer to relatively 
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recent books [2 El HI] where its different versions can be found. We start, liowever, 
with writing down necessary formulas in a form convenient for us. 

A differential equation of order N can, of course, be rewritten as a special system 
of A'^ differential equations of the first order. A consideration of first order systems 
without special assumptions on their coefficients gives more general and transparent 
results. A large part of the paper is written in terms of solutions of first order 
systems which implies the results about solutions of differential equations of an 
arbitrary order as their special cases. 

Let us briefly discuss the structure of the paper. In Sections 2 and 3 we collect 
necessary formulas for solutions of first order systems. They are used in Section 4 
for the construction of the integral kernel R{x, y, z) of R{z). In particular, we obtain 
a new representation for the integral 

' R{y,y,z)dy (1.19) 

■i 

where the points xi^X2 G R are arbitrary. Then passing to the limit xi — > — cxd, 
X2 — i- +00, we prove the trace formula (|l.lip for the coefficients Vj, j — I, . . . ,N, 
with compact supports. A construction of solutions of equation (jl.4p with asymp- 
totics (|1.13p is given in Section 5. Here we again first consider a general system of A^ 
differential equations of the first order. Finally, in Section 6 we give the definition 
of the normalized Wronskian for operators H with arbitrary short-range coefficients 
and extend the trace formula to the general case. At the end we prove that the 
normalized Wronskian coincides with the perturbation determinant. 

1.4. We note that there exists a somewhat different approach to proofs of 
formulas of type ()1.17|) . It consists of a direct calculation of determinant ()1.15|) 
whereas we proceed from a calculation of trace (II. 2p . In this way formula ()1.17|) 
was proven in [11] for the Schrodinger operator on the half-line. In [11] the Fredholm 
expansion of determinants was used. 

A general approach to a calculation of determinants Det (/ -I- was proposed 
in the book [9] by I. C. Gohberg, S. Goldberg and N. Krupnik. In this book integral 
operators K with so called semi-separable kernels were considered. It is important 
that operators K = VRo{z) fit into this class. This approach was applied to the 
Schrodinger operator in paper [8J. 

The authors thank F. Gesztesy for pointing out references [TTl [5]. 

2. Resolvent kernel 
In this section we consider an auxiliary vector problem. 

2.1. Suppose that the eigenvalues Cji j — of an A x A matrix 

Lq are distinct. We denote by pj — {pij,p2,j, ■ ■ ■ ,PN,jY (this notation means 
that the vector pj is considered as a column) eigenvectors of Lq corresponding to 
its eigenvalues Q and by p* eigenvectors of Lq corresponding to its eigenvalues 
Cj. Recall that (pj,p^) = if j ^ k (here (■,■) is the scalar product in C^). 
Normalizations of Pj and p* are inessential, but we suppose that (Pj,pp = 1- 
Then the bases pj and p*,j — l,...,N, are dual to each other. 

Assume that a,n N x N matrix V(a;) where a; € M belongs locally to and has 
compact support. We write a; « if x lies to the left of the support of V(a;) and 
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a; ^ if a; lies to the right of this set. Wc put 

L(a;) =Lo + V(a;). (2.1) 
Consider the homogeneous equation 

u'(a;) = L(a;)u(a;) (2.2) 

for the vector-vahied function u(a;) = {ui{x), . . . ,un{x)Y . For arbitrary hnearly 
independent solutions Uj{x) = {ui j(x), . . . ,un j{x)Y of this equation, we denote 

by 



U(x) 



^ ui^i{x) Ml, 2(2;) ... -ui,Ar(a;)^ 
^2,1(2;) M2,2(a;) ... U2,n{x) 



{ui(a;),U2(a;),...,UAr(a;)} (2.3) 



\UN,lix) UN,2{x) ... un,n{x)J 
the corresponding fundamental matrix. It satisfies the matrix equation 

U'{x) = L(x)U(x). (2.4) 



It follows that 



and hence 



ddetV{x)/dx =detU(x) tr (U'(a;)U"i (x)) 

= detU(.T)trL(a;) (2.5) 



det U(x2) = exp ( / tr L{y)dy) det U(, 

J Xi 



xi) (2.6) 

for arbitrary points xi and X2- Of course det U(x) ^ for all x G R. 

We always suppose that kj -.— ReCj ^0 ioi all j — 1, . . . ,N. Let n and iV — n be 
the numbers of eigenvalues Q of the matrix Lq lying in the right and left half-planes, 
respectively. The cases n = or n = where all Q lie in one of the half-planes 
are not excluded. Let Uj{x) be solutions of equation (|2.2p distinguished by the 
condition 

Uj{x) = e'^^^pj for X « if kj > and for a; » if Kj < 0. (2.7) 

We denote by and K_ the linear spaces spanned by all solutions Uj{x) such 
that Kj > and such that Kj < 0, respectively. Clearly, dim K_|_ = n and dim K_ = 
N — n. We assume that 

K+nK_ = {0}. (2.8) 



Then all nontrivial solutions of equation p.2p exponentially grow either as x -t- +00 
or as a; — > —00. In particular, equation (|2.2p does not have nontrivial solutions 
UG L2(K;C"). 

If Ui(a;), . . . , u„(a;) and u„+i(x), . . . , un{x) are arbitrary linear independent so- 
lutions from K+ and K_ respectively, then in view of condition (|2.8p all these 
solutions are linearly independent. It is now convenient to accept the following 



Definition 2.1. Suppose that n columns of matrix (|2.3|) form a basis in the linear 
space K_|_ and other N — n columns form a basis in K_. Then the fundamental 
matrix U(a:) is called admissible. 
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Observe that for the "free" case where V(a;) = 0, we can set 

Uo(x) = {pie';i%...,pjve^"^} (2.9) 

and 

Wo(a::) =detUo(a;) = det{pi, . . . , pat} exp ( trLox) (2.10) 

because tr Lq = X^^Li Cj- Note that det{pi, . . . , Pn} since all eigenvalues of the 
matrix Lq are distinct. The inverse matrix Go(x) — Uq^{x) satisfies the relation 

GS(a;) = {pte-^"^%...,p^e-^"-}. (2.11) 

2.2. Next we consider the nonhomogeneous equation 

<p'{x)^L{x)ip{x)+i{x), i{x)^{h{x),...,fN{x)y, (2.12) 

where the vector- valued function f{x) has compact support. Let us use the standard 
method of variation of arbitrary constants and set 

(p{x) =U(a;)q(x), q(x) = {qi{x), . . . ,qNix))\ 

so that 

N 

fix) ^'^qj{x)uj{x). (2.13) 

i=i 

Here U(a;) is an arbitrary admissible fundamental matrix (|2.3p . Then it follows 
from equation (|2.4p that 

q'(a;) = g(a;) where g(a;) = G{x)f{x) and G{x) = U"i(x). (2.14) 

We are looking for a solution of equation (I2.12[) decaying (exponentially) as 

a;| oo. It is convenient to accept convention (jl.Sp on the eigenvalues Q of the 
matrix Lq. Set 

/9_i_ = min ReCj, p_ ~ min |ReCj] (2-15) 

J — l,...,n j— n+l,---,A^ 

and observe that estimates 

Uj(x) = 0(e-''±l"^l), x-^Too, 

hold for j = 1, . . . ,n and the upper sign as well as for j — n + 1, . . . , N and the 
lower sign. Taking into account (|2.13p , we see that we have to solve equation (|2.14p 
for different components qj (x) of q(a;) by different formulas. Namely, we set 



Qj (x) = - gj iy)dy, j = l,...,n, 

J X 

qj{x) ^ gj{y)dy, j = n + l,...,N, 

where gj{x) are components of g(a;). This leads to the following result. 

Proposition 2.2. Let assumption (j2.8p hold, and let (|2.3p be an arbitrary admis- 
sible fundamental matrix. Then the function 

^ />oo ^ I^X 

ifix) ^ -Y,u,{x) {G{y)i{y)),dy+ J2 u,{x) (G(y)f (y)),dy (2.16) 

j=l "'^ j=n+l 

satisfies equation (|2.12p and (p{x) — 0(e~''±l'^l) as x -> ^oo. 
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Formula p.l6p can be rewritten as 

/oo 
K{x,y)i{y)dy (2.17) 
-oo 

where the matrix-valued resolvent kernel (or the Green function) R(x, y) = {Rk,i {x, y)} 
is defined by the equality 

ri N 

Rk,i{x,y) = -^Uk,j{x)gjj{y)9{y - x) + ^ Uk,j{x)gjj{y)9{x - y). (2.18) 

j=l j=n+l 

Here 9 is the Heaviside function, i.e., 9{x) = 1 for x > and 9{x) = for a; < 0, 
and gj^i are elements of the matrix G. In the matrix notation formula (|2.18p means 
that 

R(a:, y) = ~V{x)P+V-\y)e{y - x) + U{x)P^V-\y)e{x - y), (2.19) 

where the projections P± are defined in the representation = C" © i^^N-n 
the block matrices 



P+ = 



/„ 0\ 

o h I iN-r 



Expressions (j2.18p or (|2.19p do not of course depend on the choice of bases 
in the spaces K+ and K_. Indeed, if we choose other bases Ui(x), . . . , u„(x) 
and u„+i(a:), . . . , UAr(a;), then the corresponding admissible fundamental matri- 
ces U(a;) and iJ{x) are related by the formula tJ{x) = U(a;)F where the operator 
F : — ?> commutes with the projections P±. It follows that U(a;)P-|-U^^(y) = 

U(x)P±U-Hy). 

Evidently, the resolvent kernel (|2.19p is a continuous function of x and y away 
from the diagonal x = y and 

R(x, x + 0,z) = -V{x)P+V-^{x), 

Rix,x-0,z) ^V{x)P^U-\x). 

It follows that 

R(x,x-0,z) -R(x,x-f 0,z) =1, (2.20) 
where I is the N x N identity matrix. 

2.3. The results of the previous subsection admit a simple operator interpre- 
tation. Consider the space L^(R;C^) and define the operator Hq on the Sobolev 
class Hi(R;C^) by the formula 

Ho = 9I — Lq, d~d/dx. 

If u(a;) — u{x)pj where u G H^(K), then (Hou)(x) — {u'{x) — (ju{x))pj, and 
hence the operator Hq is linearly equivalent to a direct sum of the operators of 
multiplication hy — Q , G M., j = 1, . . . , N , acting in the space L^(M). It follows 
that the spectrum of the operator Hp consists of straight lines passing through all 
points —Q and parallel to the imaginary axis. In particular, the inverse operator 
Hq^ exists and is bounded. 
To define the operator 

H = 51 - Lo - V(a;), 
we need the following well known assertion (see paper ^ by M. Sh. Birman). 
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Lemma 2.3. Let T : L^(R; dx) — > L^(M; d^) he an integral operator with kernel 

(2.21) 

IfbiO = ie + iy^'^ and 



t{^,x)^b{Oe-"^v{x) 



lim 



x+1 



\v{y)fdy = 0, 



(2.22) 



then the operator T is compact. 



If the coefficients of the matrix V(a;) satisfy condition (|2.22p . then according to 
Lemma [2.31 the operator VHq^ is compact. Hence the operator H is closed on 
H^(K;C^) and by virtue of the Weyl theorem essential spectra of the operators 
H and Hq coincide. Condition ()2.8p implies that is not an eigenvalue of H so 
that the inverse operator exists and is bounded. If the matrix- valued function 
V(a;) has compact support, then according to Proposition 12 .21 the integral kernel of 
the operator is given by formula (|2.19p . 

2.4. Let the solutions Uj{x) of equation (|2.2p be distinguished by conditions 
(|2.7p . Let us give expressions for the Wronskian W(a;) :— detU(a;) in terms of 
transition matrices T± defined as follows. For j — 1, . . . ,n and a; » or j = 
n + 1, . . . , N and x « 0, we have 



N 



(2.23) 



k=l 



with some coefficients tj^k- Set 



t2,l t2,2 



^2.9 



(2.24) 



tn,7i / 



and 



tn+l,n+l in+l,n+2 
tn+2,n+l tn+2,n+2 



(2.25) 



\ tM,n + l tN,n+2 ■ ■ ■ tj^^N ) 

Consider, for example, T4.. Using expressions (|2.23p for j — l,...,n, we see 
that for a; :» matrix (|2.3p equals 

N N 

U(x) - {^^l,fcPfce^-^...,^^„,fcPfce^-^p„+le^"+l^...,PJve';""}. (2.26) 
fc=l fc=l 

Below, by the calculation of determinants of matrices, we systematically use that 
one can add to each column another column multiplied by an arbitrary number. In 
particular, we have 



W(x) = det { ^ ti^kPke'^"'' , ■■■,Y1 ^".fcPfce^"'"' Pn+ie 



,PNe 



} 



k=l 



k=l 



:detT+ Wo(a;), a; » 0, 



(2.27) 
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where the free Wronskian \Nq{x) is given by formula ()2.10p . In view of relation 
([23]), it follows that for all a; e M 

W(a;) = exp ( tr Lox - / tr V(?/)dy) det T+ det{pi, . . . , pA?}. (2.28) 

J X 

Quite similarly, using expressions (j2.23p for j = n + 1, . . . , N and x « 0, we 
obtain that 

W(a;) = dot T_ Wo(x), x « 0, (2.29) 

and 

W(x) = exp(trLoa;+ / tr V(y)d7;) det T_ det{pi, . . . , p^v}, Va; G M. (2.30) 

This leads to the following result. 

Proposition 2.4. Let\J{x) be the fundamental matrix (12.31) where Uj{x) are the 
solutions of equation (j2.2l) satisfying conditions (j2.7p . Let the transition matrices 
T± be defined by formulas ^^^^ - ([223- Then the Wronskian \N{x) = detU(a;) 
admits representations (j2.28p and (j2.30|) . 

Putting together equalities (j2.28p and (|2.30l) . we see that 

/oo 
tr V{y)dy) detT_. 
-oo 

Assumptions (j2.8p and det T± ^ are of course equivalent. 

3. Dual problem 

Some properties of admissible fundamental matrices become more transparent if 
one considers the dual problem corresponding to the matrix- valued function 

L{x) = -L*(x). 

3.1. It follows from equation ()2.4p for U(a;) that the inverse operator G(a;) — 
U^^(a;) satisfies the equation 

G'{x) ^ -G(a;)U'(a;)G(a;) = -G(a;)L(a;) (3.1) 

which yields the equation 

ij'{x) ^Z{x)ij{x) (3.2) 

for the matrix-valued function U(a;) := U*(a;)~^. Clearly, detU(a;) so that 
U(a;) is a fundamental matrix for this equation. Proposition 13.11 below shows that 
it is admissible. Set U(a;) — {■Ujj{x)}, G{x) — {gj,i{x)}. We use below that 

= 9j,i{^) = {-iy+'mi,j{x)/\N{x), (3.3) 

where mi,j{x) is the minor of the matrix U(x) which is the determinant of the 
matrix cut down from U(a;) by removing the row with index I and the column with 
index j. 

Proposition 3.1. Let \i.j{x) be arbitrary linear independent solutions of equation 
from K+ for j = 1, . . . , n and from K_ for j — n + 1, . . . , N , and let XJ{x) be 
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the corresponding admissible fundamental matrix p.Sp . Then for all I = 1, . . . ,N 

we have 

gj^iix) ^OieP-""), x^-oo, j = n + l,...,N, 
gj^i{x) =^0{e-P+'=), x~^+oo, j = ^^'"^^ 
with positive numbers p± defined in (12.15^ . 

Proof. Let us prove, for example, the first of relations (|3.4[) . Changing if necessary 
the numeration, we can suppose that j — N which is notationally convenient. For 
X <<C and some numbers Cj^k, we have 

n n N N 

fc=l k=l k=l k=l 

It follows that 

n n 

mi^^ix) = det{ ^ Cl,fcp(''e^'=^ . . . , ^ c„,fep(')e';'-% 

fc=l A:=l 
N N 

Cn+i,kP^^e<-\ . . . , ^ CN-i,kP^^^e<-^} (3.5) 
fc=i fc=i 

where 

Pk^ = bl.fe; • • ■ ,Pl-l,k,Pl+l,k, ■ • • ,PN,kY- (3.6) 

Thus the clement with index I is removed from p^ = . . . ,PN,kY so that in 

p.Sp we take the determinant of the (A^ — 1) x {N — 1) matrix. Neglecting columns 
which are repeated in determinant (13.51) , we see that it consists of terms 

det {p^'^e?-, . . . , p«e^-, pi']^^e'^'="+^% . . . , pil,e^'=«-^} (3.7) 

multiplied by some coefficients which do not depend on x. Here the indices fcn+i, . . . , k^-i 
take the values n+1, . . . , N and kp ^ kq ii p ^ q. Evidently, expression p.7p equals 

TV 

Cj exp ( trLox ~ Qx) , tr Lq = ^ Cfc, (3.8) 

k=l 

for some j = n + 1, . . . , and a number Cj which does not depend on x. In view 
of formulas (|2.10p and (|2.29p after division by W(a::) this expression is 0(|e^^J^|) as 
a; — > — oo. Hence (|3.4p for j — N follows from (|3.3p . □ 

Let iij{x) be columns of the matrix \J{x), that is, 

U(a;) = {ui{x), . . . ,ujv(a;)}. 
Then relations (|3.4p can equivalently be rewritten as 

Uj{x) = 0{eP-'"), X ^ ~oo, j^n+l,...,N, 
u.j{x) ^ 0{e^ X —i' +00, j = l,...,n. 



(3.9) 



Let us define the resolvent kernel R(a:;, y) in the same way as in subs. 2.2 with 
L(a;) replaced by L(a;). According to relations (j3.9p the first n columns of the matrix 
U(a;) exponentially decay as a; — >■ +oo and its last N — n columns exponentially 
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decay as x — > — oo. Therefore the fundamental matrix U(a::) is admissible for 
equation p.2|) , and it follows from formula (|2.19|) applied to the dual problem that 

Ti{x,y) = -V{x)P+\J-\y)e{y - x) + l5{x)P -\]-\y)e{x - y) (3.10) 

where 

U(a;) =U*(a;)-\ P± = P^. (3.11) 

In particular, we see that 

R-(a^:y) = -R-*(y,2;)- 

This relation implicitly also follows from the results of subs. 2.3. 

3.2. Let Go(a;) — {gf}(x)} be the matrix inverse to the free matrix ()2.9p . The 
next result supplements Proposition 13.11 and plays an important role in our proof 
of the trace formula (|l.lip . 

Proposition 3.2. Let solutions Mj{x) of equation (|2.2|) he distinguished by condi- 
tions (|2.7p . and let U(x) be the corresponding admissible fundamental matrix ()2.3p . 
Then for all I — 1, . . . , N elements gj,i{x) of the matrix G{x) — XJ^^{x) satisfy the 
relations 

9jAx)-9fhx)^0{e-''+n, x^+oo, j=n + l,...,N, 
with positive numbers p± defined in (I2.15P . 

Proof. Let us use the notation introduced in the proof of Proposition 13. II We shall 
again prove the first of relations (|3.12l) for j = N. According to (j2.26p for a; ^ 0, 
we have 

N N 

mi,N{x) = det { ^ ^l,fcpi'^e';'=^ . . . , ^ t^^kp'i^ e<^\p'^l,e<-+^^ , p^^l,e<--^'-^} 

k=l k=l 

(3.13) 

where pj^,''' is vector p.6p obtained from p^ by removing the component with index 
I. Neglecting columns which are repeated in {N — 1) x {N — 1) matrix p.l3p . we 
see that 

n 

mi^Nix) = det { ^ ii,fcpi'^e^-" + ^l.A^pi^^e^"^ . . . , 
fc=i 

n 

E i",.pi'^e^- + tn,NP^;,\^^^, Pi'lie^--, . . . , P^lie^"--} . 

k=l 

This determinant is the sum of the term 

n n 

det { J2 ^l,fcPi'^e^^^ ■ • ■ . E ^«.fcPi'^e^^^ pi'|ie^"-^% . . . , P^lie^"--} (3.14) 

fc=l k=l 

and of the n terms 

n n 

det {ti,A.p(^)ef«% ^2,fcpi'^e^^^ . . . , J] ^n,fcpi'^e^^^ pi'|le^"^^^ . . . , p^liC^"--} , 

fc=l k=l 

(3.15) 
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k=l k=l 

For the free case where V(a;) = 0, we have the exact equahty 

mf^lix) = det {p^'^e';^^ . . . , p«e^"% pi'|le^"+^^ . . . , p^^l,e<"-^^} . (3.16) 

Similarly to p.27p . wc find that determinant p.l4p equals 

det T+ det {pi"e^^^ . . . , pi')e^"^ pi'|ie^"+^% . . . , pj^liC^"--} . (3.17) 

By virtue of (12.27^ expression p.l7|) divided by W(a;) equals expression p.l6p 
divided by \Nq{x). 

Thus for the proof of asymptotics p.l2p it remains to estimate determinants 
p.lSp by Ce^P+^. This is similar to the proof of Proposition 13.11 It suffices 
to consider terms corresponding to different values of index k in different sums. 
Therefore, up to some factors not depending on x, determinants p.lSp consist of 
the terms 

det {p^^e';"^ p^^y^^^. . . . , p^^l/''^-^^,p%e<-+^^, . . . , p^^,l,e<-^^} (3.18) 

where the indices fci,...,fc„_i take the values l,...,n and kp ^ kg ii p ^ q. 
Evidently, (jS.lSp equals expression (13. 8p for some j — 1, . . . ,n and a number Cj 
which does not depend on x. In view of formulas (|2.10p and (|2.27l) after division 
by W(a;) this expression is 0{e~P+^) as x — > +oo. This concludes the proof of 
asymptotics p.l2p . □ 

In view of equality ()2.1ip in terms of columns u.j{x) of the matrix U(a;) relations 
p.l2p can equivalently be rewritten as (cf. p.9p ) 

' ^ - (3.19) 
Uj(x)-p*e ^J'' = 0(e''-^), X ^ -oo, j = l,...,n. 

We emphasize that the functions p*e~^J^ exponentially grow at infinity while the 
remainders in formulas p.l9p exponentially decay. 

3.3. Let us, finally, verify that the resolvent kernel defined by equality (|2.19p 
satisfies a natural estimate 

\R.{x,y)\<Ce-P^''-y^, p = min{/9+,p_}, (3.20) 

with a constanlQ C not depending on x and y. It suffices to check that 

|U(.t)P±U-1(?/)| < Ce""*!^-^!, ±{y-x)>0. (3.21) 

We can suppose that lJ(x) is defined by formula (12. 3p where the solutions Uj{x) of 
equation (|2.2I) satisfy condition (|2.7p . 
Note that 

|U(x)P±| < Ce^^il^l, Ta;>0. (3.22) 
Passing here to the dual problem and using that p± — Pi^, we see that 

|P±U-i(2/)| = |U*(y)-ip±| = |U(y)P^| fCCe-^il^l, ±y>0. (3.23) 

Let us prove estimate p.2ip . for example, for the upper sign. If cc < and y >0, 
then ([3?2T|) directly follows from ([3^221) and (|3?23| . 



"'^Here and below C are different positive constants whose precise values are of no importance. 
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Suppose next that a; < and y < 0- According to formula (|2.18|) it suffices to 
check that 

\ukA^)9M\ < Cle^'^'^-'H (3.24) 

for all k,l — 1,. ..,N and j = 1,. ..,n. Recall that \uk.j{x)\ < C\e''^^\ and that 
\g3,i{y)\ < C\e-'^^y\ by Proposition [321 This yields and hence (|32T1) for 

X < 0, y < and the upper sign. 

Similarly, we obtain that estimate (I3.2f p is true for a; > 0, y > and the lower 
sign. Using this estimate for the dual problem we see that 

|U(j;)P_U-i(?;)| < Ce-'^-l^-^l, x>0,y>Q,x>y. 

Passing here to adjoint matrices and taking into account relations p. lip , we find 
that 

\V{y)P+V-\x)\ < Ce'^+l^-^l, x>Q,y>0,x>y. 

Interchanging x and y, we get (|3.2ip for x > 0, y > and the upper sign. 
Thus we have proven the following result. 

Proposition 3.3. Let assumption (|2.8p hold. Then the resolvent kernel (|2.19p 
satisfies estimate (j3.20p . 

Estimate (|3.20p shows that formula (|2.17p obtained for functions f with compact 
support extends to aU f e L'^{R; C"). 



4. Trace formula 

4.1. Let us consider the differential operator (|l.ip acting in the space L^(R). 
Recall that the operator Hq — i^^d^ is self-adjoint in the space L^(R) on domain 
V^Hq) which is the Sobolev class H^(E). The spectrum a{Ho) of Hq coincides with 
[0, oo) for N even and with M for N odd. If the coefficients vj, j = 1, . . . ,N, of 
operator (|1.14p belong locally to L^, then the operator H = Hq + V is well defined 
by formula (|l.ip at least on the class (M) . If all functions Vj satisfy assumption 
p.22p . then according to Lemma [Z!3l the operator VRo{z), z ^ a{Ho), is compact. 
Therefore the operator H is closed on domain 'D{H) = T>{Hq), and by virtue of 
the Weyl theorem its essential spectrum (7ess{H) = (Jess{Ho)- The spectrum a{H) 
of the operator in C \ (Jess{Ha) consists of eigenvalues (not necessarily real) of 
finite multiplicity which might accumulate to aessiHo) only. 

For the construction of the integral kernel of the resolvent R(z) = {H — z)^^ of 
the operator H, we have to solve the equation 

i^'^(p^'^\x) + {V^){x) ^ z(p{x) + fix), z(^a{H), f e L\R). (4.1) 

Let us rewrite it as a system of differential equations (|2.12p . We introduce vectors 
(f = [tpi, . . . ,fNY with components ipj = ip'-^^^\ f = i^{0,0, ... ,0,/)* and set 



/ 


1 . 


. 


0\ 





. 


. 








. 


. 


1 


\i^z 


. 


. 
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Vix) 



/ 







\vi{x) 
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(4.3) 



Then equation (|4.ip is equivalent to vector equation (I2.12[) with the operator h{x, z) 
defined by equahty (|2.ip . We emphasize that it now depends on the spectral param- 



eter z. Matrix (14.21) has eigenvalues Q such that C- 
eigenvectors 



z with the corresponding 



p,(z) = (i,0,...,Cf-')* 



(4.4) 

l)/2 for Imz > and 



It is easy to see that n — N/2 if N is even and n 
n = (TV + l)/2 for Imz < if TV is odd. 

If the coefficients Vj{x) have compact supports, then all results of Sections 2 
and 3 apply. Now we have Uj{x, z) — {uij{x, z), . . . , UM.j{x, z))* where the func- 
tions Uj{x,z) := ui^j{x,z) satisfy homogeneous equation ()1.4|) and uj*"' ^\x,z) — 



Ukj{x,z) for k — 2,...,N. Hence fundamental matrix (|2.3p takes form (|1.3p . 
This matrix U(a:, z) is admissible if the functions ui{x, z), . . . , u„(x, z) belong to 
L^(R_) and the functions Un+i{x, z), . . . ,U]si{x, z) belong to L^(M-|_). The Wron- 
skian \N{x,z) — detU(a;,z) satisfies equations (|2.5p and ()2.6|1 where trL(a;,z) = 



Wat (a;). This yields ()1.9|) . Formula (|2.17p means that 



(x) 



R],N{x,y,z)J{y)dy. 



In particular, we have 

Rj,n(x, y, z) = 9^"^i?i,Ar(a;, y, z), j 



,iV. 



(4.5) 



Condition (|2.8p is equivalent to the assumption that z is not an eigenvalue of the 
operator H . Thus, Proposition 12.21 entails the following result. 

Proposition 4.1. Suppose that z ^ cr(H). Let the matrix Il{x,y, z) — {Rj^k{x,y, z)} 
be defined by formula ()2.19p where XJ{x, z) is an admissible fundamental matrix. 
Then the resolvent R{z) = [H — z)^^ of the operator H is the integral operator with 
kernel 

R{x, y, z) = i^Ri^N{x, y, z). 
According to formula (|2.18p Proposition 14. 1 1 implies that 



R{x,y,z) = -i^y^ui 



R{x,y,z) =i 



N 



N 

E 



,jix, z)gj,N{y,z), x<y, 
ui,j{x,z)gj^N{y,z), x>y, 



(4.6) 



It follows from relations p.20p and (|4.5p that, for N > 2, the function R{x, y, z) 



as well as its derivatives Ri''\x, y, z), k = 1, 
X and y while 



,N — 2, wee continuous functions of 
Ri^''^\x,x- 0,z) - Ri'^-^\x,x + 0,z) 
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The case = 1 is trivial. Although in this case the kernel R{x, y, z) is not 
a continuous function, the difference R{x,y,z) — Ro{x,y,z) is continuous and its 
diagonal values equal zero. 

4.2. Let us find a convenient expression for integral (|1.19p . Since V does not 
depend on z, differentiating equation (12.41) in z, we have 

U'(a;, z) = L(x, z)\J{x, z) + Lo(z)U(a;, z). 

Here and below the dot stands for the derivative in z. Using formula we now 

see that 

d{G{x, z)\J{x, z)) /dx ^G'{x, z)\J{x, z) + G{x, z)ij'{x, z) 

^G{x, z)to(z)lJ(x, z) A{x, z). (4.7) 
Next, we calculate the matrix A{x,z). It follows from formula (|4.2p that 



Lo(2;) = i 



N 



and hence 



B(a:, z) := Lo(z)U(a;, z) ~ i 



where (cf. (|1.3p and (12. 3p ) uij ~ uj. Since A = GB, this yields the relation 





... 


"1 







... 







Vi 


... 


0^ 






/ 





... \ 

























Wl,2 ■•• Ui^N / 



N 



'^^9j,lh,k = gj,NbN,k = i^gj.NUi^k 



(4. 



1=1 



for the matrix elements aj,k = aj,k{x, z) and hj,k — bj,k{x, z) of the matrices A and 
B. 

Putting together equalities (14. 7p and (|4.8p . we obtain the relation 

d ^ 

i^ui^k{x,z)g.j^N{x,z) = — gj,K^: z)ui^k{x, z) 
and, in particular. 



1=1 



N 



i^Uij{x,z)gj,Nix,z) = — gj^i {x, z)ui,j {x, z) 



(4.9) 



1=1 



Comparing formulas ()4.6p and (|4.9p , we get two representations for diagonal values 
of the resolvent kernel: 

, n N 



R{x, X, z) 



dx 



J2J29jAx^z)^i,j{x,z), 



3=1 1=1 
N N 



R{x,x,z) = — Ygj^i{x,z)ui^j{x,z). 



j=n+l 1 = 1 
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Integrating the first of these representations over an interval {xi,x) and the second 
over an interval (x, X2), we arrive at the following intermediary result. 



Proposition 4.2. Under the assumptions of Proposition^^ for all xi,X2,x G 
the representation holds: 



N N 



/ Riy, y,z)dy = - ^Yl i^^ z)'^i.3 i^^ z) 



N N n N 

- XI Y9jAx2,z)ui^j{x2,z) + YYgj^i{xi,z)uij{xi,z). 

j=n+l 1=1 j=l 1=1 



(4.10) 



Let us consider the first term in the right-hand side of (|4.10p . Since 
ddetlJ{x,z)/dz = detU(x,z)tr (U(x, z)"^U(x, z)) , 

we see that 

N N 

XX(7j-i(a;,z)M,,,(x,z) -tr(G(a;,2)U(x,.2)) ^ \N{x, z)-^\N{x, z) (4.11) 

3 = 1 1=1 

where the function \N{x,z) = detU(a;,z). Observe that according to (11.91) this 
expression does not depend on x which is consistent with formula (I4.10p . 

For = 2, identity (|4.10p reduces to formula (1.26) of paper [7]. For iV > 2, it 
is probably new. Identity (|4.10p allows us to considerably simplify the calculation 
of Tr {R{z) — Ro{z)) compared to the presentation of book for N = 2. This is 
essential for an arbitrary N. 

4.3. Now we suppose that the solutions Uj{x,z) of equation p.2p are distin- 
guished by condition (12.71) which yields condition (|1.6p on the solutions Uj{x, z) of 
equation (|1.4p : \N{x,z) is the Wronskian ()1.7p . Different objects corresponding to 
the "free" operator Hq — i^^d^ will be endowed with index (upper or lower). 
For the free case, we put u''^\x) = e^^^, j = 1,...,A^. Let \Jq{x,z) be the cor- 
responding fundamental matrix p. 91) with the eigenvectors Pj{z) defined by (|4.4p . 
and let \Nq{z) be its determinant (jl.Sp . The normalized Wronskian A{x,z) is de- 
fined by formula (jl.lOl) . We denote by gj,i{x,z) and gj'^^{x,z) matrix elements of 

the matrices G(a;, z) — U^^(x, z) and Go(a;, 2) = U^^(a;, z), respectively. 

For the proof of the trace formula, we combine Propositions l3.2l and l4.2l Indeed, 
let us subtract from equality (|4.10p the same equality for the resolvent Ro{z) — 
(Hq — z)~^ and consider 

{R{y,y,z)-Ro{y,y,z))dy. (4.12) 



xi 



Since Cf' — i^ z^ for alH = 1, . . . , A^, we have 

uiA^.z) = 4^]{x,z) = didf^e^n/dz = i''N-\~''+'-\l - I + x<:j)e^'^ 

if either j — 1, . . . , n and x <K orj — n + 1, . . . , N and x ::>> 0. According 
to condition (jl.Sp . it now directly follows from relations (|3.12p that, for all I = 
l,...,iV, 

g,,iix,z)uiA^,z)-g%\x,z)ufUx,z) = Oixe'^''^+''-^\^\) 
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if either j — 1, . . . , n and x — t- — oo ot j ~ ji + 1, . . . , N and x — > +oo. Therefore 
using equahty (|4.10p we see that the contribution to (|4.12p of the terms depending 
on xi disappear in the hmit xi — > — oo and the contribution of the terms depending 
on X2 disappear in the hmit X2 — >■ +oo. Thus taking into account equahty (j4.11l) . 
we obtain the following result. 

Theorem 4.3. Suppose that the coejjicients Vi,...,vn of the operator H have 
compact supports. Then for z ^ cr{H), the limit in the left-hand side exists and 

lim / {R{y,y,z) - RQ{y,y,z))dy = -A{x,zy^dA{x,z)/dz (4.13) 

where x is arbitrary. 

4.4. The left-hand side of (I4.13P can of course be identified with the trace 
of the difference R{z) — Rq{z). To that end, we first verify inclusion p.lSp . We 
proceed from the following well-known result (see, e.g., survey [4] by M. Sh. Birman 
and M. Z. Solomyak). 

Proposition 4.4. Suppose that 

/oo foa 
(l + C')"|6(0l'dC < OO, := {l + x^ f\v{x)fdx <oo (4.14) 

-oo J —oo 

for some a > 1/2. Then the integral operator T : L'^{M.;dx) — > L^(M;d^) with 
kernel ()2.2ip belongs to the trace class 6i and its trace norm satisfies the estimate 
ll^^llsi l£ CBV where the constant C depends on a > 1/2 only. 

Now it is easy to prove the following result. 

Lemma 4.5. Under assumption (jl.l2p inclusion (jl.lSp holds. 

Proof. Let us proceed from the resolvent identity 

N 

R{z) - Ro{z) ^ -Ro{z)VR{z) ^ -'^Ro{z)vjd^-^R{z), z^a{H), (4.15) 

j=i 

where Vj is the operator of multiplication by the function Vj{x). The operators 
d^R{z) are bounded because the operator HqR{z) is bounded. Proposition 14.41 
implies that Ro{z)vj G ©i ff > 1. Indeed, let $ : L^{R;dx) -> L'^{R;dC) be the 
Fourier transform. Then the operator $i?o(^)wj has integral kernel (I2.2ip where 
b{0 = (27r)-i/2(^JV _ ^)-i g^^^ ^^2;) = Vj{x). By virtue of (|n^ condition (pi^ 
is satisfied in this case. If iV = 1, then we can use that the operator \vi\^/^Ro{z) 
belongs to the Hilbert-Schmidt class. Thus all terms in the right-hand side of (|4.15p 
belong to the trace class. □ 

Since the kernel R{x,y,z) (and Ro{x,y, z)) is a continuous function of x,y and 
the limit in the left-hand side of (|4.13p exists, we see (see, e.g., [13], Proposi- 
tion 3.1.6) that 

/oo 
{R{y,y,z)-Ro{y,y,z))dy. (4.16) 
-00 

Putting together formulas (I4.13P and (|4.16p . we get the trace formula (|l.lip for the 
coefficients vi, . . . ,vn with compact supports. 
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5. Integral equations 

Here we consider differential equation (|1.4|) with arbitrary short-range coeffi- 
cients. Actually, we follow the scheme of Section 2 and first consider a more general 
equation (|2.2p . 

5.1. As usual, we suppose that the eigenvalues j — 1,. . . ,N, of an x A^ 
matrix Lq are distinct and do not lie on the imaginary axis. We set 

= j = i,...,A^, (5.1) 

where pj are eigenvectors of Lq and the vectors p* form the dual basis. We have 
P2 = Pj, Pj-Pfc = if fc, and 

TV 

i=i 

Let a matrix L(x) be given by formula ()2.1|) where we now assume that 

VeL\R±). (5.3) 

We shall show that, for all j = 1, . . . , A^, equation p.2|) has solutions u'^\x) such 
that 

uf'\x) = e'^^^(pj + o(l)), X ±oo. (5.4) 
Thus we construct solutions of (|2.2p both (exponentially) decaying and growing at 
infinity. We emphasize that our construction of the functions Uj^\x) (of Hx)) 

requires condition (j5.3p for x G M+ (for x G M_) only. Functions u^^-'(a;) will be 
defined as solutions of integral equations which we borrow from the book (see 
Problem 29 of Chapter 3). For definiteness, we consider the case x — J> — oo and put 

Set 

K,(x)= E Pme^"'"0{-x) - E P™e^""^W> (5-5) 

where = Re Cm ■ It follows from relations (|5.2p that 

K^- {x) = LqKj (x) - 5{x)l, (5.6) 
where 5{x) is the Dirac function. We also use the estimate 

|K,(x)|<C,( E e«'"-0(-x)+e'=^-0(x)), (5.7) 

m:K7Ti>Kj 

which is a direct consequence of definition (|5.5p . In particular, we see that 

|Kj(x)| < Cj-e"^"^, VxgR. (5.8) 

Let xx be the characteristic function of an interval X. Consider the integral 
equation 

/oo 
(a; - y) V(y)x(-oo,a) {y)^j {y)dy, x < a, (5.9) 
-oo 

for a function Uj(a;) = Uj(x; a) depending on the parameter a which will be chosen 
later. If then the first sum in (|5.5I) is absent. In this case we 

can omit X{-oo,a){v) in (|5.9p so that (|5.9I) becomes a Volterra integral equation. 
However (|5.9p is only a Fredholm equation for other values of j. Suppose that a 
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function Uj{x) satisfies the estimate u.j{x) — 0(6"^^) as x — >■ — oo and equation 
(|5.9p . Differentiating ()5.9p and using ()5.6|) we see tliat 



Uj(2;) = Oe'^^'=Pj+V(a;)x(-oc,a)(a;)uj(a;)-Lo / Kj(x-2/)V(y)x(-oo,a)(2/)uj(y)dy. 

— OO 

(5.f0) 

Putting together equations (I5.9P and (jS.lOp we find that a solution Uj (x) of integral 
equation (|5.9I) satisfies also the differential equation 

u'(a;) = LaUj{x) + V{x)x(~oo,a){x)uj{x), 



which reduces to equation (|2.2p for a; < a. 
Let us set 

Uj{x;a) — e''^^Wj{x;a), x < a, (5-11) 



and rewrite equation (I5.9p as 

w,(x;a) = p, - r K,(a; - y)e-<^^--y^Viy)Wjiy;a)dy. (5.12) 

— OO 

By virtue of assumption (|5.3p and estimate f|5.8p we can choose the parameter a 
such that 



/a pa 
\Kj{x-y)e~'^'^^-y^V{y)\dy <C \V{y)\dy<l, Vx G 
-OO J — OO 



(5.13) 



Then equation (|5.12p can be solved in the space L°°((— oo, a); C^) by the method 
of successive approximations. 

This result can also be reformulated in the following way. Let Qj{a) be the 
integral operator with kernel 

Q,{x,y) =K,{x - y)e-^^^--y'>V{y) (5.14) 

acting in the space L°°((— oo, a); C^). Then 

w,(a) = (/-Q,(a))-ip, (5.15) 

where the inverse operator exists because ||(Q)j(a)|| < 1. 

Clearly, the function Uj{x; a) defined by formula (I5.1ip satisfies differential equa- 
tion (|2.2[) for a; < a. Since a solution of a differential equation of first order is 
determined uniquely by its value at one point, it suffices to require equality (|5.1ip 
only for one x < a. Then the corresponding solution can be extended to all a; G M. 
Now we are in a position to give the precise definition. 

Definition 5.1. Let wj"^(-;a_) G L°°((-oo, a_); C^), j = 1,...,N, he the so- 
lution of equation ()5.12p where a = a_ is a sufficiently large negative number. 
We define \x;a^) as the solution of differential equation (|2.2p which satisfies 

condition (|5.1ip for some (and then for all) x < a. The solutions u^'*'''(a;; a+), 
j = 1, . . . , N , are defined quite similarly if (— oo, a_) is replaced by (a+, oo) where 
a+ is a sufficiently large positive number. 

It remains to verify asymptotics (|5.4p for the function Uj{x). According to 
()5.1ip and (|5.12p it suffices to check that the integral in the right-hand side of 
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()5.12|) tends to zero as x —oo. Using estimate (|5.7p and the inclusion Wj e 
L°°((— cx), a_); C^), we see that this integral is bounded by 



C r \Yiy)\dy + C ^ T, 



e('-^-^,)i^-y)\Y(y)\dy. (5.16) 



The first integral here tends to zero as a; ^> — cx3 by virtue of condition (|5.3p . Each 
of the integrals over (x, a) can be estimated by 

g(K„-.,W2 r |v(y)|dy+ r^^\V{y)\dy. 

Jx/2 Jx 

Since Km > Hj, this expression tends to zero as x — — oo by virtue again of condition 

Thus we arrive at the following result. 

Proposition 5.2. Let assumption (|5.3p hold, and let a+ (a-) 6e a sufficiently large 
positive (negative) number. Then, for all j — 1,...,N, the functions Uj^\x;a±) 
(see Definition lb.l^ satisfy equation (|2.2p and have asymptotics (|5.4p . 



Solutions of equation (|2.2p are of course not determined uniquely by asymptotics 
(|5.4[) . In particular, the solutions u^^''(a:;a±) generically depend on the choice of 
the parameter a±. 

Let Uj{x; a, r) — u"j^\x; a±,r) be the function constructed above for the cut-off 
coefficient Vrix) = X{-r.r){^)'^{x)', thus function ()5.15p is now replaced by 

MVj{a,r) = (/-Qj(a)x(-^,^))"^Pj. (5.17) 

Since 

lim ||Qj(a)(l - X{-r,r))\\L^{-oo,a) = 0, 

we see that Uj{x;a,r) — >■ Uj{x,a) as r ^ oo for all fixed x < a. This relation 
extends to all a; e M because solutions of differential equations depend continuously 
on initial data. Therefore Proposition 15.21 can be supplemented by the following 
result. 

Lemma 5.3. Under the assumptions of Proposition^^ let\i"^^\x;a±) and\i"^\x;a±,r) 
be the solutions of equations (|2.2p with V(a;) and Vr(a;), respectively, specified in 
Definition \5.1[ Then for all j = I, . . . , N the relation 

lim u'f^\x;a±,r) = u'^\x;a±) (5.18) 
holds uniformly in x on compact intervals of R. 

5.2. If a function Uj^\x) satisfies equation (|2.2p and has asymptotics (|5.4p . 

then adding to u'^\x) a solution with a more rapid decay (or less rapid growth) as 
X ±oo we obtain again a solution of equation (j2.2p with the same asymptotics 
(|5.4p . It is natural to expect that this procedure exhausts the arbitrariness in the 
definition of u^^\x). The precise result will be formulated in Lemma 15.51 
The following assertion is almost obvious. 

Lemma 5.4. Suppose that solutions u^^-*, . . . , u^-* of the differential equation (|2.2p 
have asymptotics (j5.4p as x — > ±oo. Then for each of the signs "±" the functions 
u[^', . . . , uj^'' are linearly independent. 
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e-^'(p, + o(l)) = Cj-.'^^'^CP' + 0(1)) + «(e"'') 



Proo/. It follows from that 

det{u['^\x), . . . , u^^ (x)} = det{pi, . . . , pn} exp ( ^ Cjx) (l + o(l)) 

as X — ±cxD. Since det{pi, . . . , Pa} 7^ 0, this expression is not zero for sufficiently 
large ±x. □ 

Lemma 5.5. Suppose that solutions u^^'', . . . , of the differential equation (|2.2p 
/lawe asymptotics (j5.4p as x ^ ±00. Let u^-^' 6e an arbitrary solution of (|2.2I) wif/i 
asymptotics (j5.4p as a; -H- ±cx). T/ien necessarily 

uf\x)^nf\x)+ 4?"l^^(^) (5-19) 

±(k1— Kj)<0 

/or some numbers c)^i . 

Proof. As before, we set Uj(a;) = ''(a^)- According to Lemma [5]3] we have 

A 

(=1 

with some numbers Cj^i. Therefore it follows from asymptotic relations (j5.4p that 

fj < "\^J J 

Kl<Kj 

as X — —00. Since the vectors pi, . . . , p^r are linearly independent, this relation 
implies that Cjj ~ ii I j and Cj,j — 1. Thus equality (j5.20p leads to (I5.19p . □ 

5.3. Let us return to differential equation (ll.4p with coefficients satisfying the 
assumption 

v,eL\R±), j = l,...,N, (5.21) 
only. Define as usual the matrices Lo(z) and V(x) by formulas (|4.2I) and (|4.3I) . 
Now ^j^ = i^ z and the vectors Pj{z), j = 1, . . . ,N, are given by formula (|4.4p . It 
is easy to control the dependence on z of matrices (j5.ip . 

Lemma 5.6. Elements p[,''j(z) of the matrices Pj{z), j ~ 1, . . . ,N, obey the rela- 
tion 

pgW = 0(|Cl'-'), \Cr^\z\^^. (5.22) 
Proof. Obviously, the basis dual to Pj{z) consists of the vectors 

p*(z) = (Cj,i, CJ^2C7^ • ■ • ' Cj^aCj"^^^) 

where the numbers Cj^i do not depend on |z|. It now follows from equality (14. 4p 
that p^] (z) = cj, zC,^"'C~'+' , which yields ^M- □ 

The next step is to control the dependence on z of matrices (|5.14p . 

Lemma 5.7. Elements q^j{x,y, z) of the matrices Qj{x,y,z), j — 1, . . . ,N, admit 
the estimate 

\qi'){x,y,z)\<C\C\'-''\m{y)\ (5.23) 
with a constant C not depending on x, y and z. 
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Proof. It follows from (|5.5p and ()5.22p that elements sl^]{x, z) of the matrix 'Kj{x)e~'^^'^ 
satisfy the estimate (cf. (|5.8I) ) 

|sg(x,z)|<C7|C|'=-'. 
This directly implies (|5.23p because, by definition (|4.3p . 

9m (a;, y, 2;) = -i^Sfei (a; - 2/' (y) ■ 

□ 

Let u^^^(a;, z; a±) (recall that u^^^ — (w^^-\ . . . , w^v ])*) ^''^ solution of equa- 
tion p.2p specified in Definition 15.11 Then the function 

u^^\x,z;a±) -.^ u['^^{x,z;a±) (5.24) 

satisfies also equation (|1.4I) and according to (|4.4I) asymptotics (|5.4p imply asymp- 
totics (|1.13p . Therefore Proposition 15.21 and Lemma [531 entail the following result. 

Proposition 5.8. Let assumption (|5.2ip hold, let |z| > c > and let a+ — a+(c) 
(a_ = a_(c)) &e a sufficiently large positive (negative) number. Then for every j = 
l,...,iV the function u^^\x^ z; a±) determined by Definition \5.i\ and equality (|5.24p 
satisfies equation (|1.4p and /las asymptotics (jl.l3p as x — ?> ±oo. Moreover, the 
corresponding solutions u'j^\x, z; a±,r) of equation (|1.4p wi/i cut-off coefficients 
X{-r,r){x)i^k{x), k — 1, . . . ,N, satisfy the relation 

lim u^-^''(a:, z; a±, r) = M^-^''(x, z; a±) (5.25) 

uniformly in x on compact intervals of M. This relation remains true for N — 1 
derivatives of the functions u^^\ 

By definition ()5.5p . the kernels Kj(a;, z) depend analytically on z except on the 
rays where ReC; — KeQ for some root Q ^ Q of the equation C,^ — i^ z. In 
addition to the condition z ^ a{Ho), this excludes also the half-line z < for even 
N and the line Rez = for odd N. Hence the same is true for the functions 
u'"j^\x, z; a±) if |z| > c > 0. Thus, for fixed x and a±, the functions Uj^\x,z;a±) 
are analytic functions of z if |z| > c > 0, Imz ^ for N even and if Imz 0, 
Rez ^ for odd. On the rays where Re^ — Re^j, the limits of u^j^\x, z;a±) 
from both sides exist but differ, in general, from each other by a term which decays 
faster (or grows less rapidly) than e''^^ as x ^> -l-oo or x —oo. 

5.4. Integral equations (15. 9p turn out also to be useful (even for functions Vj{x) 
with compact supports) for a study of asymptotics of the solutions m^^^ (x, z; a±) of 
differential equation (|1.4p as |z| — cxo. We choose the sign " — ", fix the parameter 
a = a_ and index j and drop them out of notation. 

Consider system ()5.12p of N equations for components Wk{x, z), k = 1, . . . , N , 
of a vector-valued function w(x, z). Set Wk{x,z) = (^''^^Wk{x, z) and take into 
account equality (|4.4p . Then we obtain for Wk{x, z) a system 

Wk{x,z) = l-^C}^'' j qk^i{x,y,z)wi{y,z)dy, x < a, (5.26) 
1=1 

where the elements qk,i of the matrix Q satisfy inequality (|5.23p . In particular, the 
operator in the right-hand side of (|5.26p is uniformly bounded as |z| — oo. 
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Assume additionally that vn{x) ~ 0. Then according to (|5.23|) the norm of the 
operator in the right-hand side of (|5.26|) is 0(|C|~^) as \z\ oo. Therefore for 
sufficiently large |C| system (|5.26p can be solved in the space oo, a); C^) by 

the method of successive approximations and Wkix, z) = l-t-0(|C|^^), k — 1, . . . ,N . 

As we have already seen in the proof of Proposition 15.21 the solution of system 
(j5.26p necessarily has the asymptotics Wk{x,z) = 1 + o(l), fc = 1, . . . , TV, as x ^> 
— oo. Define as usual u{x, z) as a solution of equation ()1.4p such that u{x, z) = 
e^^?Z'i(x, z) for x < a. Thus we obtain the following result. 

Proposition 5.9. Let assumption ()5.21|) hold, and let vn — 0. Fix arbitrary a±. 
Then for all j — 1, . . . ^ N and all sufficiently large \z\ equation (|1.4|) has solutions 
u'j^\x, z; a±) with asymptotics (|1.13p as x ^ ±oo and such that 

(a;,z;a±) = e?^-"(l + 0(|zri/~)), |z| oo, (5.27) 
for all ±{x — a±) > . 

Remark 5.10. If assumption (|5.2ip is true for both signs, then we can set a — +oo 
in equation (j5.26p . For sufficiently large |z|, such an equation can again be solved 
by the method of successive approximations. 

6. The Wronskian and the perturbation determinant 

6.1. Let us define the Wronskian W(a::) for differential equation (|2.2I) where the 
matrix-valued function Jj^x) is given by formula (j2.1l) and V(a;) satisfies assump- 
tion (|5.3p (for both signs) only. To justify the definition below, we start with the 
following observation. 

Lemma 6.1. Suppose that both sets of solutions Uj and Uj of the differential equa- 
tion (j2.2p have asymptotics (|5.4p as x ^ — oo for j — 1, . . . ,n and as x ~^ +oo for 
j = n -|- 1, . . . , N . Then for all x 

det{ui(a;), . . . , \\m{x)} = det{ui(x), . . . , \in{x)}. (6.1) 
Proof. Let us proceed from Lemma 15.51 To simplify notation, we suppose that 

Kl > . . . > K,i > Kn+l > ■ ■ ■ > KN- (6-2) 

First, we check that for alH = 1. . . . , n, 

det{ui(a;), . . . ,un{x)} = det{ui(a;), . . . , u/(a;), u/+i (x), . . . ,UAr(a::)}. (6.3) 

For I — 1 this equality is obvious because necessarily ui(a;) — ui(a;). Suppose that 
(|6.3p is true for some I. Then using ()5.19p we see that the left-hand side of ()6.3p 
equals 

det{ui(x),...,u;(a;),u;+i(a:;) -I- ^ Q+i,„u„(a;), u/+2(a;), . . . , ujv(a;)}. 

Km>Kl + l 

Since according to (j6.2p the contribution of the sum over > ki+i equals zero, 
this yields relation (j6.3p for ^ -I- 1 and hence for all Z = 1, . . . , n. 
Quite similarly, we can verify that for alH = A^, . . . , n -|- 1 

det{ui(a::), . . . , u„(a;), u„+i(a;), . . . , UAr(a;)} 

= det{ui(a;), . . . , u„(x), u„+i (x), . . . , u/(x), u/+i (x), . . .,un{x)}. (6.4) 

Putting together equalities (|6.3p and (|6.4p for I — n, we get (|6.ip . □ 
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Now we are in a position to define the Wronskian \N{x). 

Definition 6.2. Let Uj{x), j = 1, . . . ,N, be arbitrary solutions of equation (I2.2p 
with asymptotics (|5.4p as x — — oo if j = 1, . . . , n and as a; — > oo if j = n + . . . , N. 
We set 

W(x) = det{ui(x),...,UA,(a;)}. (6.5) 

Recall that solutions of equation (|2.2|) with asymptotics (|5.4|) exist according 
to Proposition 15.21 Although they are not unique, according to Lemma 16.11 the 
Wronskian \N{x) does not depend (up to a numeration of eigenvalues on a 
specific choice of such solutions. In particular, we have 

\N{x) ^det{u[^\x;a-),..., u^^'^ ix;a-), u|,+\ (a;; a+ ) , . . . , uj^^ ix;a+)} (6.6) 

where a+ (a„) are sufficiently large positive (negative) numbers and the solutions 
u'^\x; a±) are constructed in Proposition 15.21 

Of course. Definition 16.21 applies if the matrices T^oiz) and V(a;) are given by 
formulas (|4.2|) and (|4.3p . respectively. In this case the Wronskian W(x, z) depends 
analytically on the parameter z ^ a{Ho). Indeed, if additionally Imz 7^ for N 
even and Re z 7^ for odd, then this fact directly follows from the analyticity 
of the solutions u^^^\x, z; a±), j — l,...,N (see subs. 5.3). Moreover, according 
to Lemma [6.11 the Wronskian \N{x,z) is continuous (in contrast to the solutions 
Uj^\x, z;a±)) on the critical rays where Re(i = Re(j for some Q ^ Q. There- 
fore its analyticity in a required region follows from Morera's theorem. Evidently, 
W(a;, z) = if and only if z is an eigenvalue of the operator H. 

6.2. Let us return to the trace formula (|1.11|) established so far for coefficients 
Vk{x), k = 1, . . . , N , with compact supports. Suppose that assumption (|1.12p holds. 
Then condition (|5.21|) is satisfied for both signs. Let us approximate Vk{x) by the 
cut-off functions X{-r.r){x)vk{x). We denote by ' (x, z; a_, r), j — l,...,n, and 
by Wj^''(x, z; a+, r), j — n + I, . . . , N , the solutions of equation (II. 4p with the 
coefficients X{-r,r)Vk determined by Definition 15.11 and equality (|5.24|) . Let us use 
formula (|6.6|) for the Wronskian \Nr{x, z) for equation ()1.4|) with cut-off coefficients 
X{-r,r){x)i^k{x)- Then it follows from relation (|5.25p that 

lim W^(x,z) = W(a;,z). (6.7) 

In view of analyticity in z of these functions we also have 

lim VV^(a;,z) = W(x,z). (6.8) 

Set Hr = Ha + Vr where the operator Vr is defined by formula (|1.14p with the 
coefficients X(-r,r)^fe- Let us write down formula (j4.13p for the operator Hr and 
pass to the limit r 00. By virtue of (j6.7p and (|6.8p the right-hand side of (j4.13p 
converges to the corresponding expression for the operator H. It is possible to 
verify that the same is true for the left-hand side of ()4.13p . We shall not however 
dwell upon it and establish the trace formula in form (jl.lip . 

Using the resolvent identity 

N 

R{z) - R.r{z) = - ^ Rr{z){V - Vr)R{z), Rr{z) = {Hr - z)-\ 
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we see that for z ^ o{H) 

N 

\\R{z) - Rr{z)\\e, < C\\Ro{z){V ~ K)i?oWlle, < Ci \\Ro(z)v,{l - Xr)\\e,- 

According to Proposition 14.41 there is (for > 2) the estimate 

\\Ro{z)v,il-Xr)\\l,<cf \v,{x)\^{l + x^rdx, a>l/2, 

J I j; I >r 

whence 

Hm \\Rr{z) - Riz)\\e, =0. 

) — ^oo 

Thus, using trace formula for cut-off perturbations Vr and passing to the 

hmit r oo, we deduce it for V. This leads to the following result. Recall that 
the normalized Wronskian A.{x, z) is defined by formula ()1.10p . 

Theorem 6.3. Under assumption ()1.12p the trace formula (II. lip holds for all 
z a{H). 

If inclusion (|1.18l) holds, then equation (|1.16p is satisfied for a generalized per- 
turbation determinant 

5,„(z) = Det (/ + (z - zo)i?(zo)yi?o(^)) (6.9) 

where zq ^ '^(H). It is easy to sec that D^giz) is the usual perturbation determinant 
for the pair Ro{zo), R{zq) at the point (z — zo)~^. Of course, equation (|1.16p for 
a function D{z) fixes it up to a constant factor only. We note that for different 
"reference points" , generalized perturbation determinants are connected by the 
formula Dz^^z) — D za{zi)~^ D z^iz) ■ Moreover, if VRq{z) G 6i, then Dza{z) = 
D{zo)~^ D{z) where D{z) is the perturbation determinant (see formula (jl.lSp ) for 
the pair Hq, H. 

Comparing equations (|l.lip and (|1.16p we see that for all a; e M and all zq ^ (j{H) 
A(x, z) = C(xo,zo)exp ( - / VN{y)dy)Dzo{z) 

where the constant C{xo, zq) ^ does not depend on x and z. 

6.3. Suppose now that vjy ~ 0. Then the Wronskian W(a;, z) =: W(z) does not 
depend on x, and it is easy to deduce from Proposition 15.91 that 

\N{z)=\Nq{z){1 + 0{\z\-^/^)), \z\^oo. (6.10) 

For the proof, it suffices to choose a+ < 0, a_ > and use asymptotics (|5.27p 
at X = 0. As a side remark, we note that according to ()6.10p the set of complex 
eigenvalues of the operator H is bounded. 

It follows from (|6.10p that the normalized Wronskian (|1.10p satisfies the relation 

A(z) = l + 0(|z|-i/^), |z|^oo. (6.11) 

Since, by Proposition 14.41 VRo{z) G &i, the perturbation determinant is correctly 
defined by formula p. 151) and (see book [TU]) 

lim Det (/ + VRoiz)) = 1. (6.12) 

I Im z\^QO 

Comparing equations (jl.lip and (I1.16p . we obtain that 

A(z) = CDet(/ + yi?o(z)) (6.13) 
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for some constant C . Moreover, taking into account relations (|6.11|) and (|6.12p . we 
see that C = 1 in (|6J3l) . 

Let us formulate the result obtained. 

Theorem 6.4. Suppose that vn — and that the coefficients Vj , j — 1, . . . , N — 1, 
satisfy assumption ()1.12p . Then A{x, z) —: A(z) does not depend on x and relation 
(|1.17|) is true for all z ^ a{H). 

6.4. Finally, we note that for a derivation of the trace formula ()l.lip the ap- 
proximation of Vj by cut-off functions X{-r.r)Vj is not really necessary. From the 
very beginning, we could work with functions Vj satisfying assumption (|1.12|) only. 
Then the text of Sections 2, 3 and 4 remains unchanged if, for all = 1, . . . , N , the 
functions e^^^Pj are replaced for x <K by u^~-'(2;;a_) where a_ is a sufficiently 

big negative number and they are replaced for x ^ by \iY\x]a+) where a+ 
is a sufficiently big positive number. In particular, the definition of the transition 
matrices in subs. 2.4 can be given in terms of the solutions u"^^\x; a±). 

However, a preliminary consideration of coefficients Vj with compact supports 
seems to be intuitively more clear. 
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